Two Jupiter-Mass Planets Orbiting HD 154672 and HD 205739 
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D ■ ABSTRACT 

We report the detection of the first two planets from the N2K Doppler planet 
^ ! search program at the Magellan telescopes. The first planet has a mass of 

Msini = 4.96Mjup and is orbiting the G3 IV star HD154672 with an orbital 
period of 163.9 days. The second planet is orbiting the F7 V star HD205739 
! with an orbital period of 279.8 days and has a mass of Msinz = 1.37Mjup- Both 

planets are in eccentric orbits, with eccentricities e = 0.61 and e = 0.27, respec- 
tively. Both stars are metal rich and appear to be chromospherically inactive, 
^ ! based on inspection of their Ca II H and K lines. Finally, the best Keplerian 
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model fit to HD205739b shows a trend of 0.0649 m s ^ day ^, suggesting the 
presence of an additional outer body in that system. 

Subject headings: Stars: Planetary Systems — Techniques: Radial Velocities — 
stars: individual (HD 154672, HD 205739) 



Introduction 



Ongoing Doppler radial velocity surveys of nearby stars have detected over 200 extrasolar 
planets in the past decade (IButler et al.ll2006l ). Those surveys focus on late F, G, K, and 
M dwarfs within 50 pc and most of the planets they have found to date are more massive 
than Saturn, and are presumably gas giants. Recently, several Neptune-mass an d lower mass 



planets have been detected , most of them with orbital periods of a few davs (IButler et al. 



Udry et al 



2006 



2007 



Endl et all 120081 ) 



20041: iMcArthur et al.ll2004j: ISantos et al.l l2004l: iBonfils et allbood . I2OO7I : iRivera et al 



2005 



As the search for new extrasolar planets continues, Doppler surveys now look through 
a broad parameter space, including long-period Jupiter analogs, very low mass planets in 
short-period orbits, multiple planetary systems, and new planets around stars with spectral 
types tha t extend beyond th e ones traditionally been searched, i.e. KOV to F8V. The N2K 
program (jpischer et al. 2005) is a Doppler survey with distributed observing campaigns at 
the Keck, Magellan and Subaru telescopes, and is primarily aimed at increasing the number 
of known hot Jupiters. Because of their proximity to the host stars, the atmospheres of 
hot Jupiters can be as hot as 2000 K, resulting in detectable emission at IR wavelengths. 
This makes these sho rt period planets ideal targets for spaceborn follow up to observe ex- 
oplanet atmospheres (IHarrington et al.l 120071 ). especially when they transit their host stars 
and undergo a secondary eclipse so tha t small differential changes in emission from the 
star-planet system can be measured (e.g. ICharbonneau et al.ll2005l : iKnutson et al.l 120071 ). It 
is also important to have a sample of hot Jupiters that is large enough to provide mean- 
i ngful constraints on the formation, migration, and evolution mechanisms of these planets 
JPord fc Rasioll2006h . 



The N2K program searches a fresh sample of the "next 2000" stars, not included in other 
current Doppler surveys, by extending the search radius to distances of 100-110 parsecs from 
the Sun. This project is intentionally biased towards metal rich stars to exploit the correla- 



tion between formation of gas giant planets and high stellar metallicity (ISantos et al.ll2004 



Fischer fc Valentill2005l ). The search strategy is optimized for the detection of Jupiter-mass 
planets with orbital periods shorter than 14 days by obtaining radial velocity measurements 
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on three consecutive nights. Variabihty in the velocity of the host star over this short 
timescale flags the star as a candidate host for a short-period planet. However, gas giant 
planets in longer orbits are also identified and detected with additional observations. As an 
example, the N2K programs at the Keck and Subaru tel escopes have already detected seven 



plane ts with periods between 21 days and 3.13 years (IFischer et al.l 120071 : iRobinson et al. 



2003). 



Here, we report the detection of the first two planets from the Magellan N2K program. 
This program has been underway at the Magellan telescopes at Las Campanas Observatory 
in Chile since 2004, and includes ~ 300 FGK metal rich stars at distances between 50 and 
100 parsecs. The new planets are orbiting the stars HD 154672 and HD 205739. Both host 
stars showed significant radial velocity scatter in the first year of observations and were 
initially flagged as hot Jupiter candidates. However, follow-up observations over subsequent 
years revealed the presence of longer period planets and a trend on the radial velocity curve 
of HD 205739 that continues to increase after over 3.5 years. 



Characteristics of the Host Stars 



HD 154672 is classified as a G3IV st ar, with apparent magnitude V = 8.21 and color B- 
V = 0.71 (Hipparcos Catalog: lESAlll997l ). The Hipparcos parallax of the star is 15.2 ± 1.11 
milliarcseconds, placing it at a distance of 65.8 ± 4.8 parsecs. The distance and the apparent 
magnitude of the star give an absolute visual magnitude of My = 4.12. The bolometric 
luminosity of the star is L^oi = 1-88 Lq, where we have included a bolometric correction 



of -0.09 derived from the empirical transformations of IVandenBerg fc Clenj (120031 ). using 
the effective temperature, surface g ravity and metallic i ty of the star. Our high resolution 
spectroscopic analysis, described in IValenti fc Fischerl (120051 ). yields Te// = 5714 ib 45 K, 
log g = 4.25 ± 0.08, ^;sinz = 1.0 ± 0.5 km s-^ and [Fe/H] = +0.26 ± 0.04 for HD154672. 
The radius of the star derived from the Stefan-Boltzmann relation and the values of the 
luminosity and effective temperature above is 1.39 -R0. We have also der ived a stellar 
mass of 1.06 Mq, a radius of 1.27 Rq and an age of about 9.3 Gyr using the iTakeda et al. 



(120071 ) grid of evolutionary models, base d on the Yale Ste l lar Ev olution Code and tuned 
to the uniform spectroscopic analysis of IValenti fc Fischerl (120051 ). The resultant log g is 



4.26lo;g5, agreement with the results of the spectroscopic analysis. The uncertainties in 



these parameters correspond to a 95% credibility interval using Bayesian posterior probability 
distributions. The derived stellar parameters of HD154672 are summarized in the second 
column of Table 1. 



The second star, HD 205739, is a F7V with V = 8.56 and B-V = 0.546 (Hipparcos 
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Catalog; lESAl Il997l ). The Hipparcos parallax of the star is 11.07 ± 1.12 milliarcseconds, 
placing it at a distance of 90.3 ±9.1 parsecs. This sets the absolute magnitude of HD 
205739 to My = 3.78, and its bolometric luminosity to Lhoi = 2.3 Lq. The value of L^oh 
in cludes a bolometric correcti on of -0.03 derived from the same empirical transformations 
of IVandenBerg fc ClemI (120031 ) mentioned above. Our spectroscopic analysis yields a T^ff = 
6176 ± 44 K, log g = 4.21 ± 0.08, vsmi = 4.5 ± 0.5 km s-\ and [Fe/H] = +0.187 ± 0.05. 
Using the effective temperature and stellar luminosity with the Stefan-Boltzmann relation, 
we calcu late the radius of th e star to be 1.33 -R0. The stellar mass, radius and age derived 
from the iTakeda et al.l (120071 ) grid of evolutionary models are in this case 1.22 Mq, 1.33 Rq, 



and 2.9 Gyr, and log g is 4.29l[j:[j^. The parameters of HD205739 are summarized in the last 
column of Table 1. 

Finally, the Ca II H & K lines of HD 154672 and HD205739 (Figure 1) indicate that 
their chromospheric activity is low. We can therefore reject activity as the cause of the 
observed radial velocity variations of the stars. Based on these observations, we adopt a 
conserva tive upper lirnit to the expected jitter (or astrophysical noise) of the stars of about 
4 m s"^ JWright et aDl2004l ). 



3. Doppler Observations and Keplerian Fits 

Doppler observations were carried out between mid 2004 and February 2008 at the 
Magellan Clay telescope using the MIKE spectrograph (Bernstein et al. 2003), with the 
addition of an iodine cell behind the spectrograph's entrance slit to model the instrumental 
profile and to set an accurate reference wavelength scale (Butler et al. 1996). The typical 
signal-to-noise of our spectra is about 130, producing photon-limited uncertainties of 2 to 
4 m s^^. Two additional sources of noise are present in the data. The first one is the 
stellar jitter estimated in §2. The second source of noise is systematic instrumental errors, 
which for MIKE has a root mean square (rms) deviation of 5 m s~^, as derived from a 
subset of observed stars that appear to have stable radial velocities over the time span of 
the observations. A sample of stable stars measured with MIKE is presented in Figures 1 
and 2 of Minniti et al. (2008). 

We obtained a total of 16 radial velocity measurements for HD 154672 and 24 measure- 
ments for HD 205739. Those measurements are summarized in Tables 2 and 3, including the 
observation dates and the radial velocity formal uncertainties introduced by photon-limited 
noise. The data are represented in Figures 2 and 3. 

For each data set, we modeled the radial velocities to fit single planet Keplerian orbits 



- 5 - 



using a Levenberg-Marquardt fitting algorithm. In the case of HD205739 it was necessary 
to include an additional variable linear trend to best reproduce the observed radial velocity 
variations. The parameter uncertainties of each best model fit were estimated by running 
1000 Monte Carlo trials on each data set, where the model result of each trial was subtracted 
from the individual data points and the residual velocities scrambled and added back to 
the velocities predicted by the models, before running a new trial fit. The adopted final 
uncertainties of each parameter are derived from the standard deviation of all the model 
trials. 

The parameter values of the best Keplerian model fit for each target are summarized in 
Table 4. The best model for HD154672 has an orbital period of 163.94 ± 0.01 days, radial 
velocity semi-amplitude Ki = 225 ± 2 m s~^, and orbital eccentricity e = 0.61 ± 0.03. The 
rms to the fit is 4.36 m s~^, with a reduced a/x^ = 1-60 relative to the radial velocity formal 
uncertainties. Adopting a stellar mass of 1.06 Mq, we derive a planetary mass of 4.96 Mjup, 
and an average relative separation for the system of 0.597 AU. The radial velocity data are 
plotted in Figure [21 together with the best Keplerian model fit. 



In the case of HD 205739, the best model has an orbital period of 279.8 ±0.1 days, 
radial velocity semi-amplitude = 42 ± 3 m s~^, and orbital eccentricity e = 0.27 ± 
0.07. The radial velocities also show a substantial linear trend of 0.0649 ± 0.0002 m s~^ per 
day that continues after over 3.5 years of observations. The rms of the data to the fit is 
8.67 m s~^, with a reduced \fx^ = 2.13 relative to the radial velocity formal uncertainties. 
Adopting a stellar mass of 1.22 Mq, we derive a planetary mass of 1.37 Mj^p, and an average 
semimajor axis for the system of 0.896 AU. The data with the best Keplerian model fit are 
represented in Figure [31 As seen in the figure, the residuals to the fit of a single planet plus 
a long term trend still appear large, showing points that deviate about 2a from the average 
~4 m s~^ precision of the individual data points, however, an analysis of the periodogram of 
those residuals reveals no significant peaks, so we can not discard nor confirm the presence 
of additional shorter period planets in this system with the current dataset. 

The amplitude of the observed radial velocity variations for each star is 10-100 times 
larger than the uncertainties of the individual radial velocity measurements, what makes 
the possibility that the detected signals are caused by noise fluctuations very unlikely. We 
quantitatively assert this statement by pe rforming a false alar m probability (FAP) analysis 



of the data using the met hod described bvlMarcy et al.l (l2005l ) (see §5.2), with the inclusion 



of possible linear trends ([Wright et al.l 120071 ). Figure [H shows the result of 1000 FAP trial 
tests for HD205739. The FAP, i.e. the fraction of trials of scrambled velocities that yields 
lower Xa than the best reported fit, is less than 0.1%. A similar analysis of the HD154672 
data gives a negligible FAP (<C 1.0%). In this last case, the median Xa of the FAP histogram 
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is about 30, the first percently Xu is 16.1, and the minimum Xa after 1000 trial tests is 12.5. 
None of the FAP trial fits produces a Xa lower than the value reported above for the best 
Keplerian fit for HD 154672. 



Discussion 



We present two new Jovian-mass planets orbiting metal rich stars. 

HD 154672b is a fairly massive planet with a mass of Msinz = 4.96 Mj^p and a very 
pronounced orbital eccentricity of e = 0.61, that causes the planet to move from 0.23 to 
0.96 AU between periastron and apastron. The planet will therefore experience surface 
temperature changes of about 300K along its orbit, reaching a maximum temperature at 
periastron of about 600K, assuming the albedo of the planet is low. If water is present in 
the atmosphere of HD 154672b, it could transition between gaseous and liquid phases along 
the planet's orbit. 

When placed in the eccentricity vs. orbital period parameter space diagram of known 
exoplanets illustrated in Figure [5], HD 154672b shows an orbital eccentricity larger than 90% 
of the discovered planets, and is only the seventh planet found with an orbital period shorter 
than 300 days and an eccentricity larger than 0.6. Of t he other six planets, four have bee n 
found to be either in multiple-planet systems ( HD 74156: iNaef et al.ll2004l : iBean et al.ll2008l ). 



to have a brown dwarf compani on fHD 3651: iMugrauer et al. 



20061 ). to be part of a wide 



stellar binary system (HD 80606: iNaef et al.ll200ll ). or to present a large radial velocity tren d 
induced by a distant body associated to that system (HD 37605; IWittenmver et al.l 120071) . 
The high eccentricies in those cases can be explained by either Kozai oscillations (IKozaill 19621 ) 
or chaotic evolution of planetary orbits in multiple systems. Another planet in this subgroup, 
HD 17156, has been recently reported to have a large orbital axis misalignment with respect 
to the st ellar rotation axis, t hat can be best explained by gravitational interactions with other 
planets (INarita et al.l 120071 ). There is however no evidence for additional objets associated 
with the last planet in this subgroup, orbiting HD 89744. The radial velocity curve of 
HD 154672 in Figure [2] shows no trend nor significant residuals to the fit that indicate the 
presence of other massive objects in that system. 

HD 205739b has a mass of Msini =1.37 Mj^p, with an average relative separation of 
0.896 AU and an eccentricity of e = 0.27. In the eccentricity vs. orbital period diagram 
in Figure [5l the parameters of this planet do not seem atypical. For planets with orbital 
periods longer than 20 days, the mean eccentricity is 0.29, therefore HD205739b has an orbital 
eccentricity that is typical of detected planets that have not experienced tidal circularization. 
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The separation of HD205739b from its host star changes from 0.65 AU to 1.14 AU between 
periastron and apastron. The maximum surface temperature of this planet is expected to 
be of the order of 400K, and the amphtude of its surface temperature change will only be of 
about lOOK along the entire orbit. One peculiarity of the radial velocity curve of HD 205739 
is the presence of a pronounced trend of 0.0649 m s^^ per day, indicating the presence of 
an additional outer body in the system with an orbital period longer than the 3.5 year time 
span covered by our observations, and a radial velocity semi-amplitude greater than 35 m 
s~^. Finally, the residuals of our best fit are a factor of two larger than expected, what hints 
the possible presence of other bodies in this system. However, more observations are needed 
to confirm this hypothesis. 
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Table 1: Stellar Parameters 



Parameter 


HD 154672 


HD 205739 


V 


8.21 


8.56 


Mv 


4.11 


3.78 


B-V 


0.71 


0.54 


Spectral Type 


G3 IV 


F7 V 


Distance (pc) . 


65.8 


90.3 


Lboi/LQ 


1.88 


2.3 


[Fe/H] 


0.26 (0.04) 


0.19 (0.04) 


^e// 


5714 (30) 


6176 (30) 




0.54 (0.5) 


4.48 (0.5) 


logg (cgs) 


4.25 (0.08) 


4.21 (0.08) 



Mstar {Mq)'' . (0.97) 1.06 (1.17) (1.16) 1.22 (1.30) 
RstariReT ■ ■ ■ (1-18) 1-27 (1.37) (1.23) 1.33 (1.43) 
Age (Gyr)** . . . (6.92) 9.28 (11.44) (1.72) 2.84 (3.76) 



"Values derived from evolutionary models. 



Table 2: Radial Velocities for HD 154672 



JD-2453000 


RV 


O'RV 


(days) 


(ms^"'") 




189.7132 


-167.7 


2.9 


190.7083 


-169.7 


2.8 


191.7204 


-173.4 


3.3 


254.5062 


149.1 


2.6 


596.6893 


104.9 


3.1 


810.9097 


-31.1 


2.5 


872.8136 


234.8 


2.5 


1189.8750 


-111.7 


2.5 


1190.8402 


-83.3 


2.8 


1215.8605 


216.6 


2.5 


1216.7893 


217.8 


2.6 


1217.8725 


224.9 


2.8 


1277.7025 


41.6 


2.7 


1299.6210 


-25.6 


2.6 


1339.5574 


-172.5 


4.1 


1501.8960 


-168.9 


2.8 
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Table 3: Radial Velocities for HP 205739 



JD-2453000 


RV 




fdavs) 






189.8080 


-33.3 


4.0 


190.8643 


-27.7 


4.4 


191.8243 


-31.7 


4.0 


254.6094 


-77.0 


3.6 


550.8993 


-39.0 


4.3 


551.8692 


-57.0 


3.9 


655.6332 


-27.5 


3.2 


657.5999 


-10.4 


3.6 


658.6117 


-11.9 


3.4 


685.5468 


23.8 


3.3 


872.8832 


-46.8 


3.5 


982.7455 


38.4 


4.0 


988.7200 


38.7 


3.9 


1013.6776 


23.1 


3.6 


1066.5152 


-13.4 


3.7 


1067.5231 


-2.7 


3.6 


1216.9394 


19.3 


4.0 


1277.8176 


71.4 


5.6 


1299.8091 


55.8 


3.9 


1300.7959 


37.7 


3.6 


1338.7675 


18.7 


4.6 
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Table 4: Orbital Parameters 



Parameter HD 154672 HD 205739 

P (days) 163.94 (0.01) 279.8 (0.1) 

Tp (JD-h2450000.) . . 3045.3 (0.1) 3390.7 (0.7) 

ijj{deg) 265 (2) 301 (8) 

Eccentricity 0.61 (0.03) 0.27 (0.07) 

ii'i(ms-i) 225 (2) 42 (3) 

dv/dt [j&s-^days-^) 0.0000 0.0649 (0.0002) 

arei{AU) (0.580) 0.597 (0.617) (0.881) 0.896 (0.915) 

Msini (Mj) (4.61) 4.96 (5.36) (1.28) 1.37 (1.44) 

Nohs 16 24 

rms (ms-i) 4.36 8.67 

Jitter (m.s^i) 4.0 4.5 

reduced ^/xj l-()0 2.13 
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Fig. 1.— Call H lines of HD 154672 (top) and HD 205739 (middle), compared to the same 
spectral line region for the Sun (bottom) . The lack of emission in the line cores is consistent 
with low chromospheric activity levels. 
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Fig. 2. — Radial velocity measurements for HD154672 over 3.6 years. The minimum mass of 
the planet, assuming a stellar mass of 1.06 Mq, is Msinz = 4.96 Mjup. There is no evidence 
for additional planets in this system. 
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Fig. 3. — Radial velocity measurements for HD205739 over 3.6 years. The minimum mass of 
the planet, assuming a stellar mass of 1.22 Mq, is Msinz = 1.37 Mj^p. The best Keplerian 
fit shows a significant trend of 0.0649 ± 0.0002 ms^ per day, suggesting the presence of an 
additional outer body in the system. 
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Fig. 4. — Empirical evaluation of the FAP of the single-planet Keplerian fit plus a linear 
trend model for HD205739 reported in §3. The Xa of that fit is ~ 1.4 and is indicated by the 
vertical line in the plot. The vertical axis shows the number of trials that produce a given 
Xa- Less than 1 of the 1000 trial fits produce Xo- values lower than the original time series 
of the observations, indicating that the FAP of the reported fit is less than 0.1%. 
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Fig. 5. — Orbital eccentricity versus period diagram for known extrasolar planets. The size 
of the symbols scales with the planets' mass. The two open circles around the data symbols 
indicate the location of HD154672 and HD205739 in this diagram. 



